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We study the effects of interference on the quenching dynamics of a one-dimensional spin 1/2 XY 
model in the presence of a transverse field (hit)) which varies sinusoidally with time as h — ho cos u)t, 
with |t| < tf = -k/lj. We have explicitly shown that the finite values of i/ make the dynamics 
inherently dependent on the phases of probability amplitudes, which had been hitherto unseen in 
all cases of linear quenching with large initial and final times. In contrast, we also consider the 
situation where the magnetic field consists of an oscillatory as well as a linearly varying component, 
i.e., h(t) = ho cosuit + t/r, where the interference effects lose importance in the limit of large r. Our 
purpose is to estimate the defect density and the local entropy density in the final state if the system 
is initially prepared in its ground state. For a single crossing through the quantum critical point with 
h = ho cos cut, the density of defects in the final state is calculated by mapping the dynamics to an 
equivalent Landau Zener problem by linearizing near the crossing point, and is found to vary as \/Ti) 
in the limit of small lu. On the other hand, the local entropy density is found to attain a maximum 
as a function of u) near a characteristic scale u>o. Extending to the situation of multiple crossings, 
we show that the role of finite initial and final times of quenching are manifested non-trivially in the 
interference effects of certain resonance modes which solely contribute to the production of defects. 
Kink density as well as the diagonal entropy density show oscillatory dependence on the number of 
full cycles of oscillation. Finally, the inclusion of a linear term in the transverse field on top of the 
oscillatory compo nent, results to a kink density which decreases continuously with r while increases 
monotonically with u>. The entropy density also shows monotonous change with the parameters, 
increasing with r and decreasing with uj, in sharp contrast to the situations studied earlier. We do 
also propose appropriate scaling relations for the defect density in above situations and compare the 
results with the numerical results obtained by integrating the Schrodinger equations. 

PACS numbers: 



I. INTRODUCTION 

The phenomena of quantum phase transitions occur- 
ring at absolute zero temperature have attracted serious 
attention of scientists in recent years 1 ^. A quantum crit- 
ical point is associated with a diverging length scale (£) 
and a diverging time scale (£ r ) which satisfy the scaling 
forms £ ~ \d\~ u and £ r ~ £ z in the vicinity of the quan- 
tum critical point. Here d denotes the deviation from 
the critical point, and v, z are the corresponding criti- 
cal exponents. Following the possibility of experimental 
studies of quantum systems trapped in optical lattices 5 , 
there is a recent interest in theoretical studies of the re- 
lated models 4,5 . When a parameter of the Hamiltonian 
of the system is swept across a quantum critical point. 
The diverging relaxation time near the quantum critical 
point forces the system to undergo non-adiabatic evo- 
lution irrespective of the rate of change of parameters. 
If the system is in its ground state at initial time, non- 
adiabatic transitions are manifested in the occurrence of 
non-zero 'defects' (called kinks) and non-zero local en- 
tropy in the system. 

According to the Kibble- Zurek argument, if a param- 
eter of the Hamiltonian is varied as t/r, the density 
of defects (n) in the final state is expected to scale as 
n ~ T -vd/(yz+x) ^ wn ere d is the spatial dimensionality 
of the system&£&2.The above scaling form has been ver- 
ified for quantum spin systems quenched through criti- 



cal point o 10 ^ 11 ' 12 ^ 13 and also generalized to the cases of 
non-linear quenching 1 ^ when a parameter is quenched as 
h(t) ~ \t/T\\ a >sgn(t), for gapless systems^ and also for 
quantum systems with disorder— and systems coupled to 
external environment. Recently, a generalized form of 
the Kibble-Zurek scaling has been introduced which in- 
cludes a situation where the system is quenched through 
the multicritical points which shows that the general 
expression for kink density can be given as n ~ T -d/2z 2 ^ 
where z% determines the scaling of the off-diagonal term 
of the equivalent Landau-Zener problem close to the crit- 
ical point. 

In parallel to the studies on spin chains in con- 
densed matter physics, great progress has also been made 
in the realm of quantum optics in exploring the dy- 
namics of two level systems undergoing landau Zener 
tunneling due to oscillatory temporal variation of the 
parameter s 19 ' 20 ^ 1 ' 22 . The Landau Zener transition prob- 
abilities have been calculated for single crossing as well 
as for multiple crossings. Superposition of a linear field 
along with the sinusoidally varying energy levels gives 
rise to an altogether different situation, which has also 
been studied thoroughly in recent years 2 ^. 

In our present work we exploit the techniques used in 
the above papers to explore the dynamics of one dimen- 
sional spin 1/2 chain undergoing quantum phase transi- 
tions due to the application of an oscillatory or an oscilla- 
tory as well as a linearly varying magnetic field and com- 
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pare the results with the earlier findings. Investigation 
of the dynamics of one dimensional transverse XY model 
under repeated quenching of linearly varying transverse 
magnetic field has been carried out in a recent work 24 
and was shown that the defect density decreases in the 
reverse path though the entropy density increases mono- 
tonically with the number of quenching. However, the 
scaling of the defect density and the local entropy den- 
sity when the quantum critical point is crossed due to 
a sinusoidal variation of the magnetic field has not been 
attempted before. Sinusoidal quenching puts an upper 
bound to the initial and final times, which makes the 
process of coarse graining invalid in the present scenario. 
The resultant dynamics of the system becomes depen- 
dent on the phases of the probability amplitudes, leading 
to the occurance of constructive or destructive interfer- 
ences. The spin- 1/2 transverse XY chain 2 ^ discussed in 
this paper is described by the Hamiltonian 

H = - - J2 ( J *«+i + J v«+i + h <), (!) 

n 

where the cr's are Pauli spin matrices satisfying the 
usual commutation relations and the interactions and 
the transverse field are denoted by J Xl J y and h, respec- 
tively, with J x , J y and h all non-random. The interaction 
strength J x and J y are always time-independent whereas 
we shall use time-dependent transverse field in the sub- 
sequent sections. 

The Hamiltonian in Eq. (fTJ) can be exactly diago- 
nalized using the Jordan- Wigner (JW) transformation 
which maps a system of spin-l/2's to a system of spinless 
fermion o 25 ' 26 given by 

/ n-l \ 

c n = exp in ^ a j < J n > ( 2 ) 

where = (a* ± icJO/2 and the operator <t^ = 
2cJjC„ — 1 and the operators c„ satisfy the fermionic anti- 
commutation relations. In terms of JW fermions, the 
above Hamiltonian can be rewritten in Fourier space with 
a periodic boundary condition as 

H = - i [{Jx + Jy) cos k + h] (4c k + cL fe C_fe) 

fe)0 

+ i( J x - J y ) sin k (c' k c'_ k + CfcC_fc)}.(3) 

Diagonalizing the Hamiltonian in terms of the Bogoli- 
ubov fermions, we arrive at an expression for the gap in 
the excitation spectrum given b y 25 i 26 

e k = [h 2 + J 2 + J 2 +2h(J x + J y ) cosk+2J x J y cos2fc] 1/2 .(4) 

The gap given in Eq. ^ vanishes at h — ^f(J x + J y ) for 
wave vectors k — and tt respectively, signaling a quan- 
tum phase transition from a ferromagnetically ordered 
phase to a quantum paramagnetic phase known as the 
"Ising" transition 2 ^. On the other hand, the vanishing of 
gap at J x = J y for \h\ < (J x + J y ) at an ordering wave- 
vector ko — cos~ 1 (—h/2J x ) signifies a quantum phase 



transition, belonging to a different universality class from 
the Ising transitions, between two ferro-magnetically or- 
dered phases. 

The advantage of employing the JW transformation 
is that when projected to the two-dimensional subspace 
spanned by |0) and \k, — k), the Hamiltonian takes the 
form H — Hk where the reduced Hamiltonian is writ- 
ten in the form 

H k {t) = - [h + (J x + J y ) cos k] h 

h + ( J x + J y ) cos k i(J x — J y ) sin k 
—i(J x — J y ) sink —h—(J x + Jy)cosk ' 

where I2 denotes the 2x2 identity matrix. Therefore, the 
many-body problem is effectively reduced to the prob- 
lem of a two-level system with two levels denoted by the 
states |0) and \k,— k) which we refer to as diabatic ba- 
sis vectors. We shall denote the basis vectors as |0) and 
|fc, — k) as |1) and |2), respectively in this work for nota- 
tional convenience and refer the states as diabatic energy 
levels. 

The paper is organized in the following way: We have 
already discussed the model we are going to study. In 
section II, we discuss the case of oscillatory magnetic 
field but restrict our attention to the situation when the 
quantum critical points are crossed only once while in 
section III the possibility of multiple crossing is included. 
Section IV is used to discuss the dynamics with a mag- 
netic field which has both linearly varying and oscillatory 
components. In every situation results obtained through 
approximate analytical methods are contrasted with the 
numerical ones obtained by direct integration. Conclu- 
sion and summary of the results are presented in the last 
section. 



II. OSCILLATORY QUENCHING THROUGH A 
QUANTUM CRITICAL POINT: SINGLE 
CROSSING 

In this section, we shall study the spin chain driven by 
an oscillatory transverse field given by h(t) — ho cos uit 
from an initial time —tt/uj to a final time so that it 
crosses the gapless point only once during the course of 
evolution. The system initially prepared in the ground 
state |1) whereas the final ground state is |2). We shall 
evaluate the probability of the state |1) in the final state 
due to non-adiabatic evolution through the gapless point. 
As discussed above, the transverse XY chain Hamiltonian 
can be written as direct sum of 2 x 2 reduced Hamilto- 
nian matrices Hk- For an oscillatory transverse field, the 
reduced Hamiltonian gets modified to 

H k (t) = 

hg cosujt + (J x + Jy) cos k i{J x — Jy) sin k 

—i(J x — J y ) sink —ho cosuit — (J x + J y ) cosk 

For any given mode fc, the instantaneous energy gap of 
the Hamiltonian is minimum at a time io,fe such that 
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ho cos cuto,k + ( J x + Jy) cos k = where the diabatic levels 
cross each other. On the other hand, the energy gap 
vanishes for the wavevectors k — and k — n at a time 
cosiot = =f(J x + J y ) cos k/ho, respectively, signalling the 
quantum phase transition of Ising class. 

Denoting the probability amplitudes for the states |1) 
and 1 2) as Ci Ck fe(i), respectively, a general state 
vector in the reduced Hilbert space is witten as. 



Mt) = ci,fc(*)|i> + c 2 , fc (t)|2) 



(5) 



Henceforth, we set (J x + J y ) = J, and modulus of the off- 
diagonal terms = | J x — J y | sin k is denoted by A fe . Also 
the modulus of rate of change of the diagonal terms at 
time t = to,k (given by hoU> sinwto.fc) is denoted bt ak- 
Using the transformations 

CiJ(t) = C hk (t)e- 1 ^ dt ' [hoCOSUJt+J cosk] 

Q^{t) = C 2 , k {t)e~ r 5t dt ' [ - haCOSUt - J cosk] (6) 

we can rewrite the Schrodinger equation describing the 
time evolution of the probability amplitudes in the form 

.dCi y k(t) . n 2i f dt \h a cos wi+Jcosfcl 

* 77 = ^k^2,k{t)e J 

dt 

dC 2 , k (t) 

I : 

dt 



A fe Ci ifc (i)e~ 2,; F dt '[^ocos^t+cosfc] _ ^ 



It should be noted that for large values of [ho cos cut + 
Jcosfc], the phase factors on the r.h.s. of eq. (7) oscil- 
late rapidly in time. As a result the amplitudes C\^{t), 
C2,k(t), averaged over small intervals of time, remain ba- 
sically constant in time. On the other hand, close to 
t = tokt the phase factors assume stationary values, thus 
leading to non-adiabatic transition between the energy 
levels^^. 

In this section, we prepare the system in the ground 
state with initial conditions at Lot — — 7r, i.e., Ci 7r) = 
1 and C2,/t(— 71") = and the state evolves to t = so 
that the spin chain crosses the gapless quantum critical 
point only once. Using Eqs. (7), one can arrive at the dif- 
ferential equation describing the amplitude C\ k (t) given 
by 

d ' ( ' u '- 2i(hoCoaut + Jcosk)^^ + AlC ltk = (8) 



d?t 



dt 



with the probability of defect in the final state at t = 
given by pk — |Ci fc(0)| . The maximum contribution to 
the non-adiabatic transition probability comes from near 
the points where the energy gap between the instanta- 
neous levels is minimum. We can therefore linearize the 
term ho cos tot in the neighbourhood of £o,fe to get 

d 2 C dC 

1 -2i{-h u){t-to, k )smujto, k )—^+AlC ltk = 0(9) 



dH 



Eq. (9) resembles the standard Landau Zener transition 
problem^ for the linear quenching of the magnetic field 
with a variation of the field h = t/r e s where r e ff is given 
by the the rate of change of the diagonal terms of the 



Hamiltonian (1). Therefore, let us define a k = l/r e s — 
^(ei -e2)|to. fc = 2h co sin ojt 0j k = 2icy/h\- J 2 cos^. The 
non-adiabatic excitation probability 2 ^ is given as pk — 
|C lifc (0)| 2 = e~ 2 ^ , where lk = A 2 /|^( £l - e 2 )\ tQik , 
leading to 



Pk 



W h o- J '- 



(10) 



At this point, the natural question to ask is that for 
what values of the parameter ho and u>, the above re- 
lation of pk is applicable. Of course, we have used the 
non-adiabatic transition probability of the standard lin- 
ear Landau Zener problem where time t evolves from 
—00 to +00. The linearization near the crossing point 
employed above holds good only for small to. More pre- 
cisely, as discussed below the linearization approximation 
is applicable when the the time period of one cycle of the 
magnetic field (2ir/u>) is much greater than the Landau 
Zener transition time i^Thz^k) for a single crossing. The 
dimensionless Landau-Zener transition timc 2 -i2£ is de- 
fined as n LZ , k = V^T iZ , fc = |C 2 , fc (i)| 2 /^|C 2:fe (0)| 2 » 
|C 2 , fc (+oo)| 2 /(^|C 2 , fc (0)| 2 ), where k = y/a&. Using the 
above definition, we find that T^z,k ~ &-k/ot k in the adi- 
abatic limit (A k /a k >> 1) whereas in the non-adiabatic 



limit (A 2 k /a k « 1), T LZ . k is given as T LZ ,h ~ 1/V"fe- 
It should be noted that for the linear quenching of the 
magnetic field h(t) — t/r, a k — 1/r. 

Generalizing to the case of p eriodic quenching, 
Tlzm ~ Afe/(2oj^/^-o J 2 co&2 k), in the adiabatic 
limit while in the non-adiabatic limit, T]^z,k ~ 

1 / ( \j 2iOyJTi^ — J 2 cos 2 k) . The transitions are localized 
around to, k as compared to the time period for one cycle 
of the magnetic field if T^z,k is less than the time pe- 
riod of one oscillation, i.e., TLz,k << n/to. This means 
that for h\ >> J 2 , in the adiabatic limit, A k << h , 
and in the non-adiabatic limit, lu << ho- This leads 
to the conclusion that the equation (10) describing the 
non-adiabatic transition probability is valid only for large 
ho and small 10. In the defect density for small lu (i.e. 
lu < 7r( J x — Jy) 2 1 \/h\ — J 2 ) only the modes close k ~ 0, ir 
contribute, resulting to a kink density at t — given by 



^uj^—J^ 
tt\J x — J„\ 



\LU < 



k{J x 



Jy) 2 , 



(11) 



The analytical results for p k and hence the density of 
defects (obtained from Eq. 11) match exactly with the 
transition probabilities obtained by numerical integration 
of the Schrodinger Eqns for ho » A k ,Lu as shown in 
Figs. 1 and 2. From Eq. 11, we find that in the limit 
ho » J, the defect density shows a scaling form n ~ 
(/low) 1 / 2 which can be generalized using the Kibble-Zurek 
argument that assumes that the non-adiabatic transition 
is only dominant at a time when the characteristic time 
scale of the system is of the order of the rate of change of 
the Hamiltonian£i£i2. In the limit of large ho and small 
lu, we can generalize the above prescription to derive a 
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scaling form for the defect density for a single crossing 
the quantum critical point due to the periodic driving of 
the transverse field given by n ~ (/iow) l/d ^ 1/z+1 - ) where 
v and z are the exponents associated with the quantum 
critical point and d is the spatial dimension of the system. 




FIG. 1: pk vs k as obtained numerically (solid line) and an- 
alytically (dashed line) for ho = 20, id = 0.0003, j J x — J v \ = 
0.05, and J = 1. 
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defined at higher temperatures. The diagonal entropy be- 
comes identical to the previously defined Von Neumann 
entropy [s VN ), given by s VN = - f£ tr(p k lnp k )dk/n, 
when the off-diagonal terms in the density matrix, eval- 
uated at the final time, go to zero upon coarse grain- 
ing over k spac o 11 ! 12 ' 24 . We have checked the varia- 
tion of the diagonal entropy density, evaluated at the 
final time, with lo by numerically integrating Sd(k) = 
Pklnpk + (1 — Pk)ln(l — Pk) over all k, with p k obtained 
from eqn. (10). It is seen that the entropy attains a 
maximum near u> = loq — ir(J x — J y ) 2 j h^ln2 where the 
non-adiabatic transition probability (see Eq. 10) for the 
mode k = ir/2 is half, and falls off on either side of loo. 
It should be noted that loq closely resembles the char- 
acteristic time scale to appearing in the case of linear 
qucnchingii. 
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FIG. 3: Variation of diagonal entropy density with lo for one 
half cycle with ho — 20, | J x — J y \ = 0.05 and J = 1, as ob- 
tained by numerical integration of Sd(fc) using the analytical 
expressions of pk (eq. 10). The entropy for one half cycle 
attains maxima near lo ~ loo = it(J x — J y ) 2 /ho In 2. 
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FIG. 2: nvsu for ho = 20, \ J X - J v \ = 0.05 and J = 1. Solid 
line is found by numerically integrating p k over k, and the 
dashed line is the plot of eqn. (11). Analytical and numerical 
results match exactly for lower values of u. 



A. Entropy 

Following a recent paper by R. Barankov and A. 
polkovnikov~, we define the diagonal entropy Sd(k) for 
each k mode by 

Sd(k) = pu,klnpil,k (12) 

I 

where pu^ — (l\pk\l), Pk being the instantaneous density 
matrix of the system for the mode k. One advantage of 
using the diagonal entropy is that, it follows the thermo- 
dynamical relations as expected to be followed by entropy 



III. OSCILLATORY QUENCHING THROUGH A 
QUANTUM CRITICAL POINT: MULTIPLE 
CROSSINGS 

Let us now focus on the repeated quenching case when 
the spin chain is periodically driven through the quan- 
tum critical point. In the present section, it will be shown 
that interference plays a major role in deciding the be- 
haviour of the system, and for some choices of param- 
eters, the phases will add up destructively to make the 
tunneling probability almost zero. In order to be able to 
treat the successive Landau Zener transitions , as inde- 
pendent events, the time between two successive cross- 
ings has to be greater than the Landau Zener transition 
time for a single crossing as mentioned in the previous 
section. To attain this limit we shall once again restrict 
our study to large values of ho and small lo. The system 
is prepared in the state |1) at time t — 0. The diago- 
nal terms of the Hamiltonian for each modes k vanish, 
and consequently the gap becomes minimum, when the 



5 



magnetic field crosses the — Jcosk line as shown in Fig. 
4. 

When the system approaches the crossing points of the 
diabatic levels, the energy gap is minimum leading to 
large relaxation time and the system fails to evolve adi- 
abatically and the non-adiabatic transitions take place. 
On the other hand, away from the crossing points, the 
system follows adiabatic dynamics. Consequently, the 
evolution matrices associated with the system are differ- 
ent for close to and away from the crossing point o 21 ' 22 . 
Between the crossings, the system evolves following the 




FIG. 4: Variation of energy levels due to the application of 
an oscillatory transverse field. The gap becomes minimum 
at the points where the magnetic field (shown by oscillatory 
solid line) becomes equal to — Jcosfc. 



matrix 



e 4 ^ 



where j denotes the direction in which the system goes 
through the crossing points. The LZ crossing in the non- 
adiabatic region can be approximately described by the 
evolution matrix 



G 



LZ,j 



COS ^ 



X p i6LZ,j 



where the angle x IS given by 
sin 2 — = 1 — exp(— 27T7) 



(13) 



where 7 is defined in the previous section and j once 
again defines the direction of quenching (with respect to 
the crossing point). We have suppressed the notation k 
denoting the wave vector for the time being. Also 



1LZA 



v stoke 



'LZ,2 ~ 77 — "stoke 



(14) 



(15) 



where T(x) is the gamma function and s t o kes — ► tt/4 
or 9 stokes — > 0, as 7 -> or 7 ^ 00 respectively. If 
the system is repeatedly quenched with the sinusoidal 
field a series of Landau Zener crossings take place with 
the evolution of the system described by the successive 
application of matrices defined above. More specifically 
for one full cycle, i.e., tot going from to 27r, one can 
write the complete evolution matrix as a product of Gj 
and Glzj, given by 



— L Z, 2<J 2<~r LZ, 1<J1; 



(17) 



which can be generalized for N complete cycles to the 
form 



Gn = {Glz^GiGlz^Gx 



(18) 



The probability amplitude of the states \i) C^jv (i = F 2) 
at the final time u>t — 2Nn therefore obeys the relation 



G2.N 



{GlZ,2G2GlZ,\G\) N 



Ci,o 

C2,0 



where C\fi = 1 and C2,o 
little bit of algebra yields 



(19) 

at initial time t = 0. A 



Glz,2G2Glz,iG\ 

where 
9u 



9ii 921 
-921 9*ii 



(20) 



cos 2 ^ e -i(Si+e 2 ) _ sin 2 X e i(-e LZA +e LZ , 2 -0 1 +e 2 ) 



g 21 = S i n I CO s |( e <(fe*.i+«i-fl2) + e *(8LZ,2+e 1 + e 2 



(21) 



Denoting the probability that for mode k the system is 
in state |2) after the nth crossing by Q n ,k, we get, 



i,fe = (1 ~Pk) 



(22) 



as seen in the previous section and for one complete full 
period of oscillation 



Q2,k = 4p fc (l - p k ) sin {Ostokes + 62) 

For small anisotropy i.e., (J x 
ujy/fiQ — J 2 cos 2 k we have 9 s tokes - 
and 92 are given by 



(23) 



Jy? 
7r/4, and 



<< 



yj (ho cos ~uxb + J COS fc) 2 



2ho + Jir cos k 



(24) 



n = 



2 / y (ho cosuit + Jcosk) 2 
-2hg + Jit cos k s 



(- 



-)• 



(25) 



Subsituting Eqs. 24 and 25 in Eq. 23, we ge" 



19,20 



htokes = -| + argT(l - ij) + "/(Inj - 1) (16) 



, _ . . 9 ,2h + Jtt cos k tt 
Q 2 ,k = 4p fc (l - Pk) sin 2 (— + -). 

QJ 4 



(26) 
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<?2.fc as obtained numerically and analytically as given in 
Eq. 26 are plotted as a function of k in fig. (5) and (6). 
The numerical plot is fairly in agreement with the analyt- 
ical results. It is seen that Q^,k oscillates, with the tun- 
neling probability going to zero for many fc's, showing the 
signatures of constructive and destructive interferences. 

It is clear that for very small o>, Qi,k oscillates rapidly 
with k due to the presence of the sinusoidal term in eq. 
(26). As a result the coarse grained or average Q2M (de- 
noted by Q2k)i obtained by integrating each Q2,k over 
a small range around that k followed by normalization, 
gives 



Q 0.0008 
2,k 



Qi,k = %PkO--Pk), 



(27) 



as obtained previously for repetition under linear 
quenching 24 . It has been shown in recent work o 11 ' 12 ? 24 
that, for linear quenching, we can evaluate the transi- 
tion probabilities at the end of each cycle by using the 
coarse grained density matrix only, thereby simplifying 
the problem greatly by neglecting the off-diagonal terms 
in the matrix. Analogously, in the present case also, the 
characterstic time scale associated with the rate of change 
of the off-diagonal terms in the density matrix for a mode 
k sets the critical value of u> below which we can safely 
denote the tunneling probability by the coarse grained 
expression of Qik only, thereby yielding eq. (27)> 19 . One 
concludes that in the limit of very small u>, the time in- 
terval between two successive crossings is large enough 
to destroy the phase information in the coarse grained 
probabilities. 
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Jy\ = 



FIG. 5: Variation of Q2,k with k for ho = 20, \J X 
0.005, u> — 0.01 and J = 1. The widely spaced dashed line 
is analytical, and the numerical data points are shown on the 
closely spaced dashed line. 




FIG. 6: Variation of Q 2 ,k with k for ho = 20, | J x - J y \ = 0.005, 
lj = 0.01 and J — 1 obtained by averaging out the oscillations 
of the data shown in fig. 5. The dashed line is analytical, and 
the solid line is numerical. The smooth dotted line is the plot 
of course grained excitation probability Q 2 ,k as obtained from 
eq. (27). 
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FIG. 7: Variation of actual and course grained kink density 
and entropy density with u for ho = 20, \ J X — J v \ = 0.005 
and J = 1 for one complete cycle, obtained by numerically 
integrating Q 2 ,k and Q 2 k in eq. (26) and eq. (27) respectively. 
The actual and course grained plots for kink density match 
reasonably well, whereas we see a significant difference in case 
of entropy density, eventhough their qualitative behaviours 
show similarity for a wide range of u). The peaks of the plots 
occur near uiq. 



We shall call the diagonal matrix as U\ and the other 
as U2 which successively operate on the column matrix 
(Ci.fe, C2.fe). The angles ( and 9 are given as 



sin 



c 



N , . 2 X 

— « 4 sin — cos ( 



'LZ,l 



'LZ,2 



2) 



(29) 



To generalize to the case of many complete periods, it A 

is useful to recast eqn. (20) in the form as2i ^ = * an ~~q (30) 

Gs, , • r ' ( cos q sin a6 ^ \ where 

V-sm|e * cos I ) A = cos 2 - sin (6»i + 6> 2 ) 

e -i0/2 \ 2 

e ^/2 ) ( 28 ) +sin 2 |sin(0 LZ , 1 -0 LZ . 2 +0 1 -0 2 ) 
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and 

B = cos 2 — cos ( 
2 

X 



h + 62) 



- sin 2 - cos (9 LZ ,i - 6lz,2 



'LZ.l 



9 



LZ,2 



(31) 



The dynamics described by Eq. (28) can be understood 
by exploring the properties of the rotation matrices U\ 
and C7 2 . The role of U% is to rotate a vector about an axis 
in the x — y plane by an angle of £, whereas the matrix 
Ui brings about a rotation of the vector by an angle 9 
around the z axis only-2 2 -. If 9 is a multiple of 27r, which 
we call the resonance condition, the z-axis rotation does 
not affect the dynamics and the small oscillations of £ 
add up constructively to produce full oscillations between 
the diabatic states |1) and |2). On the other hand, if 9 
differs from a multiple of 2ir by more than (, then the 
rotations about angle £ do not add up constructively, and 
the oscillations will be suppressed, thus resulting in an 
effective rotation about an axis almost parallel to the z 
axis onlySi (see figure 9). In the present context, the 
resonance condition is given by, 

9&2(9 1 +9 2 ) = 4?rJcosfc =2 mr, (32) 



2 J cos k 



(33) 



Therefore for the resonance conditions, we can write 



C2,N 



± 



cos 



C 



L P -i4> 



sin 

cos i 



N r 



Ci,o 

C-2,0 



From this formalism it is clearly seen that for the res- 
onance conditions, after N complete cycles, we get N 
successive rotations by the small angle £. This causes 
oscillations in the probabilities of the two states with fre- 
quency given by 



n = 



2-k/uj 



2ir 



COS ( V s toke 



Pk 

2ho + Jir cos k s 



(34) 



Since C depends on the wave vector k, the proba- 
bilities for each resonant mode oscillates with its own 
characteristic frequency. As a result the kink density 
as well as the entropy density obtained by integrating 
over all modes shows an oscillatory behaviour (see fig- 
ure 10). The oscillatory behaviour of the entropy den- 
sity observed here is an artifact of retaining the phase 
information of the off-diagonal terms of the density ma- 
trix. It can be shown that in absence of phase informa- 
tion, \Ci t k(t)\ 2 , |C2,fc(F)| 2 — > 1/2 after each crossing, and 
as a result the entropy density of the system increases 
monotonically— with the number of crossings. 




FIG. 8: Graph showing the behaviour of the excitation 
probability as a function of number of full cycles for ho = 
20, (J x - J y ) = 0.005, J = 1 and w = 0.01. The dashed 
line is obtained analytically for k — 81.9521° with integral 
2 J cos k/ui = 28, while the solid line is the analytical graph 
for k = 25.8419° with integral 2J coak/to = 180. Numeri- 
cal data points shown on the dashed line corresponding to 
k = 81.9521° match exactly with the analytical results. 




FIG. 9: Graph showing the behaviour of the excitation prob- 
ability as a function of number of full cycles, obtained analyt- 
ically for h = 20, 1 J x - J v \ = 0.005, J = 1, u = 0.01, k = 80° 
and non-integral 2 J cos k/uj. As expected, the excitation 
probability varies randomly and does not differ appreciably 
from its initial value. 



IV. QUENCHING BY A MAGNETIC FIELD 
VARYING BOTH LINEARLY AND 
PERIODICALLY 

In this section, we study the defect generation in a 
transverse XY spin chain driven by a time-dependent 
magnetic field h(t) which consists of a linear part as well 
as an oscillatory part given by h(t) — t/r + h$ cos cot 
where r denotes the rate of the linear part of the quench- 
ing. In the limit ho — ► 0, the dynamics reduces to the well 
studied Kibble-Zurek problem of linear quenching while 
in the other limit of r — * 00, one should recover the re- 
sults presented in earlier sections. The presence of both 
the linear and periodic terms non-trivially modifies the 
density of defect in the final state as shown below. The 
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FIG. 10: Graph showing the behaviour of the kink density 
(solid line) and entropy density (dashed line) as a function of 
number of full cycles for ho = 20, \ J X — J y \ = 0.005, J = 1, 
and uj = 0.01, as obtained by numerically integrating the 
excitation probabilities obtained by using eqn. (19). 



reduced Hamiltonian in the present situation is given by 
H k (t) = 

e i{J x — J y ) sin k 

-i{J x — J y ) sin k —e 

with e = t/r + ho cos ut + J cos k and we shall once again 
recall the parametes ctk and as defined before. 




FIG. 11: Graph showing the behaviour of the diabatic energy 
levels with time when a linearly varying magnetic field is ap- 
plied in addition to the oscillatory term. The inclined dotted 
line is the plot of t/10, the solid line is ^ + cost, and the 
dashed line parallel to x axis is the constant — J cos k = 1. 

For a given wavevector k, the instantaneous energy gap 
is minimum at times to k such that 



+ ho cos ujto k + J cos k = 

T 



(35) 



Since coswio.fc < 1 always, so from the above equation 
we can conclude that all the io.fc's occur in the time in- 
terval — t(J + ho) < t ,fc < T {J + h a )- Also, since time 
between two successive io,fc' s is of the order of tt/lo, so we 
can estimate the number of times that the gap goes to 



minimum for any k is ~ , with the minimum number 
of times being 1. The situation is depicted in Fig. 11. 
In the adiabatic limit (A 2 /^ >> 1) the Landau Zener 



transition tim o 29 ' 33 (tlz) is given as tlz ~ ^i/ a k, and 
in the non-adiabatic case, we have tlz ~ l/y/c*k- It 
should be noted in the present case, the rate of change of 
diagonal terms ak = |jj2(- + h cosujt + Jcosk)\ tok — 
2|- — hoOJ sino;£o,fc|- Hence for our theory to be valid, i.e., 
to get widely separated non-overlapping LZ transitions, 



we need r << 



for the adiabatic and Jr « — 



for the non-adiabatic situations. 

We prepare the system in its ground state at t — > — oo 
with \Ci t k{— oo)| 2 = 1 and the probability of defect for 
the mode k in the final state at t — > +oo is given by the 
probability |Ci,fe(+oo)| 2 . For the linear as well as peri- 
odic driving, Eq. 8 when linearized around the crossing 
point t = to } k gets modified to 



d 2 d 



+ 2i(( h uj sin ujto,k)(t 

a z t t 



*o,fc)) 



dt 



0, 



(36) 



which leads to the non-adiabatic transition probability 

Pk = e i-"o-=in-t , fc (- 37 ) 
In the limit of small r and large u), we can expand the 



ff(Jx — Jy) sin k 



Fur- 



excitation probability as p k r» 1 - xpr_ hoOJ sinoJto k ■ 
ther, for — » hocusmcoto^k, we can write the expression 
for kink density as 



[1 



7r( J x — Jy ) 2 sin 2 fe 
1/t — /iqw sin wto.fc 



(38) 



The approximate equation given in Eq. (38) matches per- 
fectly with the numerical integration results (see Fig. 13). 
On the other hand, in the limit of large r and small w, 
only the modes close to k = or 7t contribute to the de- 
fect density, and by considering only the and k modes 
in £o,fc) we can arrive at an approximate analytical ex- 
pression given by 



1 



p k dk 



T y i — /io w sin (x^o.i 

2 u/a: J?v 



ZlQW Sin W^o.tt 



^\Jx Jy | 



(39) 



The kink density as a function of r for the non-adiabatic 
and adiabatic cases, as obtained from equations (38) and 
(39) respectively, are plotted in Figs. 13 and 14 together 
with the corresponding numerically obtained values. As 
expected, in case of non-adiabatic evolution, we get exact 
matching between analytical and numerical results only 
for low values of r for which - >> how sinto,fc, whereas 
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in case of adiabatic evolution, we see good agreement 
between the analytical and numerical results only when 
- is not close to ho sinuio.o or sinojto^, since around 
these values of r, the effects of h smojt .k for k ^ 0,n 
become important. 

It should be noted that for the single crossing case, 
using the similar line of arguments given in section II, 
one can propose a generalized Kibble- Zurek scaling form 
of the defect density in the final state given as 



1 



ii/d/(i/z+l)i 



n <~ ao[| houjsinujto o| 

T 

+ a 7I [\--h cjsinu;t , 7T \^/^ z + 1 '>} 

T 

where qq and are two constants. 



(40) 




20 40 60 80 100 120 140 160 180 

k 

FIG. 12: pk vs k (in degrees) for the case when gap becomes 
minimum only once, with t = 2, ho = 1, ui = 0.1, \J X — J y \ = 
0.05 and J = 10. The dashed line is analytical and the nu- 
merical results shown as data points coinside exactly with the 
analytical values. 




FIG. 13: Kink density (n) vs r for ho = 1, ui = 0.1, | J x — J y \ = 
0.05 and J = 10. The solid line is the plot of eq. (38), and 
the dashed line is obtained by numerical integration of the 
analytical expression of pk as given in eq. (37). We get exact 
matching between the two results for low r only, as expected 
from the theory. Only single crossing occurs for the range of 
r shown in the figure. 




20000 30000 40000 



50000 60000 70000 



FIG. 14: Kink density n vs r for ho = l,w = 0.000001, A = 
0.05 and J = 10. Only single crossing occurs for the range of 
r shown in the figure. The dashed line is found by numerical 
integration of eq. (37) , and the solid line is obtained by using 
eq. (39). We get exact matching between the two results 
except when i is very close to how sin uto,o or hoco sin tot ^. 



Now we concentrate on the situation of multiple cross- 
ings of the energy minima. For multiple crossings to oc- 
cur for a given wave vector k, Fig. 11 implies that there 
should exist a t = t such that - — /i wsinu;? = 0, i.e., 
sin cut = —, 



-. This is possible only if | — -^i— - 1 < 1. The 
Schrodinger equations for the probability amplitudes can 
be put in the form 



dCi.k 
dt 



Akel (4+ 2 - h ^+^cosk) C2k 



dC 2 



dt 



(41) 



(42) 



Using the relation e ±lzsinu>t = T,T=-oo J r {z)e ±lrut , 
where J r (z) is the Bessel function of first kind of order 
r, given by 



-iy 



Jm=0 m!r(m + r + l) v 2 



(?) 2m+ \ 



(43) 



we recast the equation to the form 



^^ = A k V j r { 2 J^y^+ 2Jtcosk+r ^c 2 , k . (44) 
dt — ' uj 

r— — oo 

The terms on the R.H.S. of eqn. (44) being rapidly 

varying in time, dC J t ' k attains a non-negligible value only 
when the phase is stationary. Using 

t 2 

h 2 Jt cos k + rut 

T 

1 2 Jt cos k + tot 2 (2 Jt cos k + ruor) 2 

= -:(*+ 2 ] a~t ( 5) 



we find that 



dd,, 
dt 



is non-negligible only close to t 



_ 2 Jr cos fc+ru;^ wit h r = 0,±1,±2, .... The above relation 
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implies the existence of an effective phase transition for 
each value of r. Choosing r — I (say) and denoting Ci ik 
by Ci k,i (i = 1,2), we get 



dCi 



l.k 



dt 



A k Jie 



(2.7 T COS fc+l^T) 



(46) 



Invoking upon the transformation to a new variable 
s = t + Iujt/2, we get 



dCxr 



1 -I 2 u 2 r 2 +4jT 2 Ii 

' : A k J,e- 1 — 



-iJ 2 T cos 2 /Cg^(s+ Jr cos A;) 2 £i ( 



2,1, k 



(47) 



Let us compare with purely linearly quenching case (/i = 
0), when the above equation gets modified to 



l.k 



ds 



A „— i J 2 t cos 2 k „ — is+Jr cos k) 2 fi 



(48) 



The role of periodic modulation on top of the lin- 
ear quenching is to renormalize A^, to A k with A^ — 

■ 2 up' t- 2 +4 Jt- 2 iaj cos k 

AfcJ;e _I & . Note that A k also vanishes at 

the quantum critical point for the modes k = and 7r. 

The dividing of the probability amplitudes for a given 
wave vector to different / values by using the Bessel's 
function can be visualized in the following way: The two 
energy levels for the wave vector k are assumed to consist 
of a number of sublevels 2 ^, with probability amplitudes 
of the Zth level being denoted by C\,i,k and C2,i.k- Each 
sublevel undergoes a level crossing only once through the 
course of dynamics, and for the It h tran sition for the 
mode fc, the incoming state (given by Ci-i^k) and the out- 
going state (given by Cj k) are connected by the transfer 



M, 



i.k 



D, 



n,k 

,2, 2 2 , 



f3i, k e~ 



D, 



l.k 



where D ltk = ^/pTJ' and $l,k = sgnJi(rf)^/l - pi ik e l * ! fc , 



in which pi_ 



-tttAV,^)) 2 



r\ = 2ho/uj and 4>i. k is the 



Stokes phase given by 

7r/4 + argT(l - iSi rk ) + 8i, k (ln8^ k 



n,k 



1) 



(49) 



where Si ik — tA 2 (J;(?7)) 2 /2, in terms of the gamma func- 
tion r(z). 

It should be noted that J;(f?) — * for large I, and 
^T^-ooJfi 1 !) — 1- So the transition is confined to a 
finite region, and the infinite series of recursive relation 
for I — ► oo converges to a finite value. In case of I for 
which Ji (j]) w 0, Mi >k gets reduced to an identity matrix. 
Hence taking J;(ry) = VZ > If we can write the state 

vector Ci.k= {C 1 .i 7k ,C 2 ,i, k ) as 

C;,fc = Mi tk C;-i,k 

= M lik M l _ lik ...M , k -M-i f+1 , k M- l/ , k C in ,k (50) 



where Cin,k denotes the initial condition. The proba- 
bility of excitation at infinite time is given by p k (oo) = 

ICijcMI 2 . 

Pk(oo) = |Ci, fc (oo)| 2 

= |[1 0]M lf<k M lf _ 1<k ... 

M 0tk ...M-i f+ i tk M-i ftk C in ,k | 2 (51) 

We are interested in evaluating the defect density in 
the final state in the limit of large r and hence the off- 
diagonal terms in the matrix Mi >k vanish upon coarse- 
graining, i.e., upon integration over all k. This approx- 
imation leads to the final result for the coarse grained 
non-adiabatic transition probability given by 



|Ci,i,k| 2 =pi,fc|Ci,i_i lfc | 2 + (1 -pi,fc)|C 2j i-i,fc| 5 



(52) 



after neglecting the cross terms. The effects of the critical 
points are manifested by A k vanishing at k = 0, ir and 
those modes do not evolve from their initial state. 

Comparison between kink density obtained by numer- 
ical integration of Schrodinger equation and by using the 
approximate analytical eqns. (51) and (52) for different r 
has been shown in figure (15). We also plot the diagonal 
entropy density in figures. (16) and (17). Although, the 
defect density in the final state decays with increasing r, 
the entropy density is found to increase and ultimately 
saturates as r increases. This is in sharp contrast to the 
cases of linear quenchin g 12 i 24 where it attains a maxima 
at a characteristic scale t — tq, and falls off on either 
side of To. This behaviour can be explained in the fol- 
lowing way: for very small r the system fails to evolve 
appreciably and therefore remains very close to its initial 
ordered state at the final time. For larger values of r, the 
probabilities change which introduces disorder in the fi- 
nal state leading to higher value of entropy density . For 
very large r, the linear term becomes insignificant com- 
pared to the oscillatory term in the Hamiltonian. Conse- 
quently, only the oscillatory quenching term contributes 
to the dynamics of the system, keeping the entropy of the 
system constant for a fixed value of u>. It is to be noted 
that the interaction between the levels depend on the 
value of Ji(2ho/uj) which saturates in the limit of large 
ui if ho is held fixed irrespective of the values of I and 
k. As a result, the kink density and also entropy density 
saturate in the limit of large lo, as shown in fig. (17) 



V. CONCLUSION 

In conclusion, we have studied the effects of interfer- 
ence in the quenching dynamics of a one-dimensional 
transverse XY spin chains in the presence of a time- 
dependent magnetic field h(t) = hocosujt or hit) = 
t/r + hocosuit. The system is initially prepared in its 
ground state and we estimate the defect density and en- 
tropy density in the final state following the quench using 
both approximate analytical and direct numerical inte- 
gration techniques. In all the situations, the analytical 
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FIG. 15: Graph showing the variation of the kink density with 
r for multiple crossings occurring with ho = 1, (J x — Jy) — 
0.05, (J x + Jy) = 10, and ui = 0.1 . Results obtained by 
numerically solving Schrodinger equation (dotted line) match 
reasonably well with the analytical ones (dashed line). 



FIG. 17: Graph showing the variation of the kink density 
(solid line) and entropy density (dashed line) with u> as ob- 
tained analytically for ho = 1, (Jx — Jy) = 0.05, (J x + Jy) = 
10, and t = 50. 



1 




0.2 / 

°0 100 200 300 400 500 600 

1 

FIG. 16: Graph showing the variation of the kink density 
(solid line) and entropy density (dashed line) with r as ob- 
tained analytically for ho = 1, (Jx — Jy) = 0.05, (J x + J y ) = 
10, and uj = 0.1. 



and numerical results are found to be in good agreement. 
Our observations presented in the paper are summarized 
as follows. 

Firstly, we have studied the defect density in the fi- 
nal state following a single crossing by linearizing the 
oscillatory magnetic field round the time at which the 
instantaneous energy gap is minimmum. We show that 
in the limit of large ho and small u>, the defect density 
scales as ^/uJ. The observation is supported by numeri- 
cal solution of the Schrodinger equation in the limit of 
small uj. On the other hand, the diagonal entropy density 
shows a maximum at a characteristic frequency scale ujq 
as defined in the text. Effects of interference are invisible 
in the case of a single crossing only. We do also suggest 
an equivalent Kibble-Zurek scaling relation for the defect 
density in the above situation. 

In the next section we generalize to the multiple cross- 
ing situation where the interference of the probability 
densities play a dominant role. We use two different 



transfer matrices valid close to and away from the cross- 
ing points. We show that for a full cycle of oscillation 
the results obtained for repeated linear quenching 24 when 
the off-diagonal terms of the density matrix are coarse 
grained, leading to loss of phase information which gives 
rise to constructive and destructive interferences, is a 
valid approximation in the limit of small uj. For mul- 
tiple crossings, we show that there exist resonance wave 
vectors for which the non-adiabatic transition probability 
oscillates between zero and one with the number of cross- 
ings following a characteristic /c-dependent frequency. As 
a result the defect density also exhibits an oscillatory be- 
haviour. The entropy density also shows a similar de- 
pendence on the number of crossings, which is in stark 
contrast with the linear quenching case, in which exclu- 
sion of the interference effects in the excitation probabili- 
ties causes the entropy density to increase monotonically 
with the number of crossings. It may be noted that a 
similar oscillatory behaviour is observed for the central 
spin of quantum Heisenberg chairs^ 4 .. 

Lastly, we study the quenching of the spin chain in 
the presence of a magnetic field which is varying linearly 
with time as t/r and also modulated by a periodically 
varying part hocosuit. For the single crossing case, we 
once again use the linearization method which predicts a 
defect density that is in fair agreement with the numeri- 
cally obtained result. For multiple crossings, we again in- 
voke the transfer matrix approach to evaluate the cross- 
grained defect density. In this case it has been shown 
that for large values r we can safely neglect the phase in- 
formation, and hence the effects of interference, by coarse 
graining the density matrix. Our analytical and numer- 
ical results show that the defect density decreases with 
increasing r for a given uj whereas when w is varied with 
r fixed, the defect density saturates for higher values of 
uj. The entropy density also exhibits a monotonic in- 
crease as a function of r with fixed w, an observation 
that is in sharp contrast with the linear quenching case 
where the entropy density attains a maximum at a char- 
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acteristic time scale Ttr^ . This may be an artifact of 
the integrability of the model which gets decoupled into 
independent two-level systems^. 
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